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ABSTRACT 

GRB 130925A is one of the recent additions to the growing family of ultra-long GRBs (T90> 

1000 s). While the X-ray emission of ultra-long GRBs have been studied extensively in the past, 
no comprehensive radio dataset has been obtained so far. We report here the early discovery of an 
unusual radio afterglow associated with the ultra-long GRB 130925A. The radio emission peaks at low- 
frequencies 7 GHz) at early times, only 2.2 days after the burst occurred. More notably, the radio 
spectrum at frequencies above 10 GHz exhibits a rather steep cut-off, compared to other long GRB 
radio afterglows. This cut-off can be explained if the emitting electrons are either mono-energetic or 
originate from a rather steep, dN/dE oc power-law energy distribution. An alternative electron 
acceleration mechanism may be required to produce such an electron energy distribution Furthermore, 
the radio spectrum exhibits a secondary underlying and slowly varying component. This may hint 
that the radio emission we observed is comprised of emission from both a reverse and a forward shock. 

We discuss our results in comparison with previous works that studied the unusual X-ray spectrum 
of this event and discuss the implications of our findings on progenitor scenarios. 


1. INTRODUCTION 

Gamma-ray bursts (GRBs) are the most energetic ex¬ 
plosions known in the Universe. These events exhibit 
prompt gamma-ray emission followed by what is referred 
to as afterglow emission in a wide range of wavelengths 
from X-rays to radio. Currently, GRBs are classified 
based on the duration of their prompt Gamma-ray emis¬ 
sion. Short (duration < 2 s) and long (duration > 2 s) 
bursts are believed to be a result of different progenitor 
systems (See reviews by Woosley & Bloom 2006; Nakar 
2007; Berger 2014). 

The common wisdom suggests that long GRBs origi¬ 
nate from the core-collapse of massive stars (e.g., Mac- 
Fadyen & Woosley 1999). Short GRBs, on the other 
hand, presumably arise from coalescence of two neutron 
stars (e.g., Paczynski 1986; Eichler et al. 1989). In these 
scenarios, the afterglow is explained by the fireball model 
(Piran 1999; Sari et al. 1999) where the broadband emis¬ 
sion originates from electrons in the circumstellar (or 
interstellar) medium (GSM or ISM) which are acceler¬ 
ated by a forward shock driven by the relativistic ejecta. 
While many studies (see above) focus on understanding 
their origin, GRBs also serve as natural laboratories to 
study physical processes in extreme conditions, such as 
relativistic particle acceleration. In turn, understanding 
these processes and the conditions in which they occur 
can provide clues as to the true nature of GRBs. This 
paper is related only to GRBs with long prompt emis¬ 
sion and thus we will discuss only this type of GRBs, 
hereafter. 

Recently, Levan et al. (2014) have suggested that an- 
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other class of high-energy transients may exist, with pos¬ 
sibly a different progenitor system. They pointed out 
that a handful of GRBs exhibit very long prompt gamma- 
ray emission (> 1000 s). These events are usually fol¬ 
lowed by a late-time X-ray afterglow which shows flaring 
activity as late as 10'^ s after the initial burst. Unlike 
other long GRB afterglows, the flaring activity in these 
ultra-long GRBs has especially high flux and longevity. 
Levan et al. suggest that these ultra-long GRBs may 
be a result of either a tidal disruption event (TDE) of a 
star by a black hole or a core-collapse of an extremely 
extended star (see also Gendre et al. 2013). 

On 2013, September 25, another ultra-long GRB, 
namely GRB 130925A, was discovered (Lien et al. 2013). 
Both the X-ray and radio afterglows of this event show 
unique features. In this paper, we discuss the unusual 
radio emission of GRB 130925A and discuss the implica¬ 
tions of its unique properties. We first summarize the 
details known so far from recent studies of this GRB in 
§ H In § [3] we describe our radio observations and the 
data reduction. The data analysis and modeling is per¬ 
formed in § m We then discuss our findings in § [5] and 
summarize in § |6l 

2. GRB 130925A - DISCOVERY AND RECENT 
STUDIES 

GRB 130925A was discovered by the Burst Alert Tele¬ 
scope (BAT; Barthehny et al. 2005) onboard the Swift 
satellite (Gehrels et al. 2004). Golenetskii et al. (2013) 
reported that the prompt Gamma ray emission, observed 
by the Konus-wind satellite, lasted for ^ 4500 s and had a 
fluence of 5.0 ±0.1 x lO^'^erg cm“^. Adopting a redshift 
oi z = 0.347, as measured by Vreeswijk et al. (2013), 
results in isotropic energy of Eiso ~ 1-5 x 10®^ erg. Com¬ 
pared to other ultra-long GRBs (Levan et al. 2014), 
GRB 130925A is the nearest event discovered so far and 
the most energetic (see Table [T]) . The prompt Gamma- 
ray emission was followed by a spectacular X-ray emis¬ 
sion with rapid flaring at the first 10^ seconds of the 
event. Once the flaring ceased, the X-ray emission de- 
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cayed as a smooth power-law, typical of normal long 
GRBs at this stage. However, additional uncharacteris¬ 
tic properties of the X-ray emission have been revealed. 


TABLE 1 

GRB 130925A in comparison to previously discovered 
ULTRA-LONG GRBS 


Name 

Z 

Eiso 

[erg] 

Tgo 

[sec] 

GRB 101225A 

0.85 

1.2 X 10®^ 

> 7000 

GRB 111209A 

0.67 

5.2 X 10®^ 

- 10, 000 

GRB 101225A 

1.77 

7 X 10®^ 

~ 6000 

GRB 130925A 

0.35 

1.5 X 10®® 

- 4500 


Notes - The properties of the previously discovered GRBs are 
adopted from Levan et al. (2014). 


Bellm et al. (2014) analysed the X-ray data obtained 
by the Swift, NuSTAR, and the Chandra satellites. Sur¬ 
prisingly, they found that the X-ray spectrum in the en¬ 
ergy range 1 — 20 keV, can not be modeled with a single 
power-law, as in essentially most normal long GRB af¬ 
terglows to date (however, see Starling et al. 2012, that 
report a thermal X-ray component in some GRBs asso¬ 
ciated with SNe) . Bellm et al. obtained satisfactory 
hts to the observed X-ray spectrum with several models 
including: 1. two power-law components; 2. a power-law 
-|- an absorption line at 6keV; 3. a power-law and a 
blackbody component. However, the physical interpre¬ 
tation of the models was inconclusive. Recently, Piro et 
al. (2014) analyzed additional X-ray data from the XMM 
space observatory. According to their analysis, the X-ray 
emission is dominated by a blackbody emission and that 
only a small contribution of the X-ray emission is due 
to non-thermal synchrotron emission from a traditional 
afterglow (an afterglow from an external forward shock). 
A different explanation for the late-time X-ray emission 
was suggested by Evans et al. (2014). They argue that 
the emission is a result of reflection from dust which re¬ 
sides far away from the GRB. Both Piro et al. (2014) and 
Evans et al. (2014) find that the lack of (or weak) X-ray 
emission from the external shock suggests an extremely 
low circumburst density (n < 0.1 cm“^). 

In the optical regime, Greiner et al. (2014) reported 
the detection of a flare, 300 — 400 s after the prompt 
emission. At late-times, Tanvir et al. (2013) observed 
GRB 130925A with the Hubble Space Telescope (HST). 
They detected an optical afterglow with an offset of 0"12 
from the host galaxy center. This offset may disfavor 
a TDE scenario for this event. However, Tanvir et al. 
also note that the host galaxy is disrupted and that this 
may be a sign of a recent merger. In this case, it is 
possible that a massive black hole can be offset from the 
optical light center, thus still leaving a TDE as a plausible 
scenario. 



V [GHz 


Fig. 1.— Radio spectra of GRB 130925A at various VLA obser¬ 
vation epochs. The initial spectrum shows a peculiar cut off at 
frequency > 10 GHz. Moreover, the spectral evolution suggest the 
existence of an underlying (and slowly variable) constant-flux com¬ 
ponent. However, more than a year after its discovery, the radio 
emission from GRB 130925A faded away below our detection limit 
(See Table [ 2 J 1 . 


3. VLA OBSERVATIONS OF GRB 130925A 

We observed GRB 130925A with the Karl G. Jan¬ 
sky Very Large Array (VLA) under a Director Discre¬ 
tionary Time (DDT) program (14A-435; PI Horesh). 
Our observations consist of multi-epoch observations 
starting 2.2 days after the bursfl Each observation was 
performed using a varying set of average frequencies: 
3.4 GHz (S-band), 6.1GHz (G-band), 9 GHz (X-band), 
14.75GHz (Ku-band), and 22 GHz (K-band). The vari¬ 
ous observing epochs (starting on 2013 Sep 27 UT; see 
Table [2]) were performed with the VLA being either in 
the A or B configuration. 

In all of our VLA observations, we used 3G48 as a flux 
calibrator and J0240-2309 as a phase calibrator. The 
data were then reduced using both AIPS and GASA (Mc- 
Mullin et al. 2007) standard routines. To estimate the 
accuracy of our flux calibration and to check for any cali¬ 
bration errors we compare the measured flux of our phase 
calibrator at the various observing epochs. The wide¬ 
band spectra of the phase calibrator at different times 
are consistent within the following wide-band flux cal¬ 
ibration errors: 6.3%, 1.3%, 2.5%, 5%, and 5% in the 
S-, G-, X-, Ku, and K-bands, respectively. Note that 
these calibration errors were calculated when using the 
full bandwidth in each band (2-8 GHz bandwidth) and 
have been adjusted for sub-band measurements. The full 
set of our measurements is presented in Table [5] 

4. RADIO SPEGTRUM ANALYSIS 

Figure [T] shows the observed wide-band radio spectra 
at different epochs. The spectrum from the first epoch, 
2.2 days after discovery, shows that the emission is al¬ 
ready peaking at low frequencies (~ 7 GHz) and that it is 
strikingly cut off at > 10 GHz. Both these properties are 

® Unfortunately, due to the U.S. government shutdown, the VLA 
operations ceased at a critical time. Our second epoch was con¬ 
ducted on the last night before the shutdown, but was limited to 
only low frequencies. Further observations resumed much later on, 
43 days after the discovery of the GRB. 












3 


TABLE 2 

Summary of JVLA radio observations of GRB 130925A 


Time 

days 

Frequency 

[GHz] 

Flux 

[fiJy] 

Flux r.m.s 

[liJy] 

2 

4.8 

237 

17 

2 

7.4 

298 

12 

2 

9.5 

293 

14 

2 

13.5 

146 

20 

2 

16.0 

89 

22 

2 

22.0 

104 

13 

9 

4.8 

216 

14 

9 

7.4 

214 

9 

9 

8.5 

ISO 

10 

9 

9.5 

203 

10 

43 

3.0 

133 

27 

43 

4.8 

105 

13 

43 

7.4 

87 

10 

43 

13.5 

60 

9 

43 

16.0 

62 

9 

43 

22.0 

83 

8 

58 

3.0 

99 

27 

58 

4.8 

77 

15 

58 

7.4 

68 

15 

58 

8.5 

65 

14 

58 

9.5 

56 

15 

58 

14.75 

46 

11 

74 

3.0 

83 

6 

74 

6.1 

71 

6 

74 

9.0 

56 

5 

74 

14.75 

42 

8 

94 

3.4 

59 

13 

94 

6.1 

47 

7 

94 

9.0 

34 

7 

94 

14.75 

33 

7 

94 

22.0 

33 

8 

108 

3.5 

39 

12 

108 

6.1 

50 

6 

108 

9.0 

41 

6 

108 

14.75 

44 

6 

108 

22.0 

38 

10 

499 

6.1 

< 24 

- 

499 

9.0 

< 24.3 

- 


Notes - Time is given in days in days since the burst. The 
errors presented in the table represent the rms error from 
each image. When given, limits are 3(t detection limits. Sys¬ 
tematic nux calibration errors were added according to §2. 
Both the rms error and the calibration should be combined 
in quadrature. 

not typical of normal GRB radio afterglows (see Chan¬ 
dra & Frail 2012 for a review of GRB radio afterglow 
properties), in particular the high-frequency cut off. 

Before discussing the implication of the high-frequency 
emission cut off, we test whether this cut off can be 
due to some modulation of the intrinsic flux via ex¬ 
treme interstellar scattering and scintillation (ISS). Ac¬ 
cording to Cordes & Lazio (2002), the Galactic scatter¬ 
ing measure (SM) towards the position of the GRB is 
logS'M ft! —3.69 [kpc The transition frequency 

from strong to weak scintillation is at ft 8 GHz. At higher 
frequencies only small ISS flux modulations are expected. 
Thus it is unlikely that the observed cut off in the radio 
emission above 10 GHz, is due to temporal strong scin¬ 
tillation (i.e., expected modulation < 10%). Moreover, 
even at frequencies below the transition frequency (i.e., 


the strong scattering regime), we do not observe strong 
scintillation (see S 14.3.31) . 

In the context of the fireball afterglow model (see § [1]) , 
GRB afterglow spectra are usually well described by a 
broken power-law (Sari et al. 1998). This is a result 
of the radio emitting electrons being accelerated into 
a power-law energy distribution, dN/dE oc E~p. The 
spectral behavior depends on some characteristic fre¬ 
quencies, which define the transition between the dif¬ 
ferent spectral slopes (see a detailed description by Sari 
et al. 1998). In short, at lower frequencies, the emission 
will be optically thick. At frequencies where the emis¬ 
sion becomes optically thin, the specific flux will rise as 
up to some maximum value after which the specific 
flux will decay as . At even higher frequencies, 

beyond some cooling frequency (which usually occurs at 
or above the optical regime, at early times), the specific 
flux has a steeper decline of 

The typical average observed value of the electron en¬ 
ergy power-law distribution is p ft 2 (e.g., Panaitescu & 
Kumar 2001), but with a relatively wide distribution of 
(Tp ft 0.5 (Shen, Kumar & Robinson 2006). Thus, typi¬ 
cally, the spectral slope of the optically thin radio emis¬ 
sion, above the peak frequency but below the cooling 
frequency, is rather shallow (iz“°-®). The sharp spectral 
cut off we observe in the case of GRB 130925A, raises 
the possibility that the power-law energy distribution of 
the emitting electrons is much steeper than in any other 
observed GRB to date. Another possibility is that an 
alternative model is required, such as a mono-energetic 
energy distribution. We next test both models by per¬ 
forming a minimum fitting to the data. 

4.1. The power-law energy distribution model 

Here we fit the initial observed radio spectrum with 
the common power-law afterglow model, consisting of an 
optically thick synchrotron self-absorbed emission 
and an optically thin emission which has a power-law 
spectral shape Since, as seen in Figure[Tl the radio 

emission appears to be decaying into a flat spectrum, it 
is possible that there is an additional slowly varying flat 
spectrum emitting component which we treat as constant 
over the time scale of our observations. Therefore, we 
perform the model fitting both with and without a second 
constant flux component. 

Figure [2] shows the best fit results. The fit without 
a constant flux component has a reduced i.e. 
per degrees of freedom (d.o.f), of xl = 4.2 (d.o.f = 3). 
Adding a second component of constant flux to the fit 
results in a = 6.2 (d.o.f = 2). The best fit spectral 
indexes in the former and latter fits are s = 1.4±0.1 and 
s = 1.6±0.2, respectively. Assuming that s = (p—l)/2, 
suggests that the electron energy distribution power-law 
indexes are p ft 3.8 and p ft 4.2, respectively. These are 
rather steep energy distributions, not observed so far in 
GRB radio afterglows. 

4.2. A mono-energetic synchrotron emission 

Motivated by the sharp cut off observed in the initial 
radio spectrum, we next consider an alternative model 

^ The assumption here is that the cooling frequency is above the 
r adio bands, an assumption which is supported by the data (see 

§ED 
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Fig. 2.— Best fit of the power-law model to the observed initial 
radio spectrum. The two fits are with (left panel) and without 
(right panel) an additional constant flux component (27/rJy). Both 
fits suggest a steep electron energy distribution, not previously 
observed in GRBs (see text for details). 


in which the observed synchrotron emission originates 
from relativistic electrons with a mono-energetic energy 
distribution. Since this model is rarely used for GRB 
afterglows (however see Waxman & Loeb 1999), we next 
describe it in more detail. 

The synchrotron emission from a single electron is 




^ f (^) 

rUeC^ \VcJ 


( 1 ) 


in the observer frame the emission is 


p 






(4) 


where T is the bulk Lorentz factor of the shockwave. The 
frequency Vsynch in the observer frame will be the same 
as above multiplied by a factor of T. 

In the case where internal synchrotron self-absorption 
(SSA) is dominant, the observed specific luminosity in 
the shockwave frame will be 


Lssa{v) 


47ri?3 p /1 - \ 

V V '^SSA ) 


(5) 


Here, we assumed a planar absorption and the optical 
depth is defined as 


TSSA = 0'SSA — , ( 6 ) 

77 

where assA is the absorption coefficient defined as (Ry- 
bicki & Lightman 1986): 


ass A = 


c 

STTlTleV^ 


N{E) d 
A2 dE 


[E^P,] dE, 


(7) 


where N(E) is the volumetric density of electrons with 
energy E = ^errieC^- 

The shape of a SSA spectrum therefore depends also 
on the energy distribution of the electrons, N{E). In the 
mono-energetic case, N{E) = No6{E — Eq) (also assum¬ 
ing constant spatial density). In this case, equation 7 is 
reduced to: 


ass A = 


S'Kiy'^ 


No 


dP, 

dE 


Eq 


2P.(Ao)jVo 

Eo 


( 8 ) 


The specific flux in this case can therefore be easily cal¬ 
culated using the following properties: B, 7 e, Ue, and 

R. 

Finally, the observed flux can be reduced to the simple 
form (see Waxman & Loeb 1999): 


where 


C = 


^synch 

1 - 


1 — [dlnF{x)/dlnx] 
A = A-K-iemeiR/Dsf, 


(9) 

( 10 ) 

( 11 ) 


in the shockwave frame, where c is the speed of light, B 
is the magnetic field strength, and me are the electron 
charge and mass, respectively, and 

pCO 

E{x)=x 74 ( 5 / 3(0 dC, (2) 

J X 

where K is the modified Bessel function. The syn¬ 
chrotron frequency, Vsynch is defined as 


^synch 


^lldeB 
AirmeC ’ 


(3) 


where 7 e is the Lorenz factor of the electrons (in the 
shockwave frame). In the relativistic case, the emission 
measured by an observer will be beamed and therefore 


and Ds i s th e luminosity distance. 

As in § 14.11 we performed the model fitting twice, with 
and without an additional constant flux component. The 
fitting of these two cases resulted in xl = 1-7 (d.o.f = 2) 
and xl = 6.2 (d.o.f = 3), respectively, where the best 
fits are presented in Figure [31 From a x^ point of view, 
the mono-energetic model with a secondary constant flux 
component is slightly preferable (although not signifi¬ 
cantly) to the power-law model from § 14.11 

4.3. Properties of the emitting medium 

We derive the electron density, the magnetic field 
strength and the radius of the emitting regions for two 
cases: 1. An isotropic sub relativistic emitting sphere. 2. 
Emission from a relativistic jet. 
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Fig. 3. — Best fit of the mono-energetic model to the observed 
initial radio spectrum. The two fits are with (left panel) and with¬ 
out (right panel) an additional constant flux component (95/rJy). 
The mono-energetic -|- a constant flux component provide the best 
fit to the data. 


4.3.1. The sub relativistic isotropic case 

Here we assume that by the time of the observation 
the bulk motion is sub relativistic, and that the emis¬ 
sion is isotropic. In addition we also assume some ratio 
between the energy density of the electrons {Ee) to the 
energy density of the magnetic field {Eb)- Adopting the 
preferred model of mono-energetic electrons and assum¬ 
ing equipartition {E^ = Eb), the best fit to the data 
at the first observing epoch requires a radius of Riso ~ 
8.5 X 10^® cm, a large electron density of Ue ~ 400 cm“® 
with 7 e ~ 45, and a magnetic field strength of B ~ 0.61 
Gauss. The minimum required energy in this scenario is 
quite large, Emin ~ 2.3 x 10"^® erg. The fact that we ob¬ 
serve only the optically thin synchrotron emission from 
the third epoch and on does not allow us to obtain esti¬ 
mates of the above model parameters in these additional 
epochs. 

4.3.2. A Relativistic Jet 

Recently, Barniol Duran, Nakar & Piran (2013), pre¬ 
sented a straightforward method for calculating the 
above properties if the emission is originating from a rel¬ 


ativistic jel0. Adopting their prescriptions and assum¬ 
ing equipartition we derive the values of the emitting re¬ 
gion properties using our measurements of the peak flux 
and frequency at our initial epoch of observation. We 
find that, at that time, the jet is only mildly relativistic 
with a bulk Lorentz factor of T Ri 4. The Jet radius is 
Rjet ~ 1-3 X 10^^ the electron density (assuming con¬ 
stant density) rig ~ 11 cm“®, the magnetic field strength 
is B ~ 0.05 G, the electron Lorenz factor is je ~ 245, 
and the minimum kinetic energy is Ek ~ 1.7 x lO'^^erg. 
The relatively low value of T should not be surprising 
given that the initial Lorentz factor is estimated to be 
To 20 — 37 (Greiner et al. 2013). 

4.3.3. A constraint on source size via scintillation 

As was already discussed above, radio flux variations, 
especially at low frequencies, can be observed due to ISS. 
The ISS flux modulations depends on the angular size of 
the emitting source. In turn, the detection or lack of such 
modulations can be used to constrain the source size. In 
our initial observing epoch, we observed the source at 
4.8 GHz twice within ^ 2 hours. We did not detect any 
strong variation in the source flux between these two ob¬ 
servations. Thus we assume that the ISS modulations 
are within our measurement errors. Conservatively, we 
estimate then that the ISS modulation is < 20%. Us¬ 
ing the equations from Walker (1998) for refractive and 
diffractive scintillation, we find that the source angular 
size should be > 3 x 10^^cm (radius > 1.5 x 10^^cm), 
comparable to the source size we derived above. 

5. DISCUSSION 
5.1. Electron energy distribution 

In §[4]we found that the energy distribution of the radio 
emitting electrons can be fitted with either a very steep 
power-law or a mono-energetic energy distribution. None 
of these distributions have been observed in GRB radio 
afterglows to date. Previously observed GRB afterglows 
have exhibited a power-law energy distribution with a 
power-law index of p « 2. Theoretical studies (see Piran 
2004 and references therein) suggest that electrons are 
accelerated at the GRB shock front via the Fermi pro¬ 
cess into a power-law distribution with a power-law in¬ 
dex ofp ~ 2.2 —2.3, in agreement with past observations. 
A steep spectrum is possible if the synchrotron cooling 
frequency is below the observed frequency. However, the 
best fit of each of the two models in ? l4.3l results in a cool¬ 
ing frequency of i/g > 10^^ Hz, well above the observed 
radio frequencies. Thus producing the observed steep 
energy distribution may require an alternative particle 
acceleration mechanism. Another recent challenge for 
current particle acceleration models in GRBs has been 
made by Wiersema et al. (2014) who detected circular 
polarization in the optical afterglow of GRB 121024A. 

Mono-energetic electrons have been observed in other 
astrophysical sources, such as the galactic center (Lesch 
& Reich 1992; Duschl & Lesch 1994). In solar flares, 
the initial particle acceleration is believed to be a result 
of magnetic reconnection which creates a mono-energetic 

® The steepness of the electron energy distribution has only a 
small effect on the values of the derived properties (Barnio Duran; 
private communication) 
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soft X-rays emission (see Benz 2008 for a review). More¬ 
over, some solar flares also exhibit steep optically thin 
radio spectra (e.g. Nita, Gary, & Lee 2004). Magnetic re¬ 
connection in accretion disks has also been suggested, as 
a possible explanations for intra-day variability in AGNs 
and Blazars (e.g., Lesch & Phol 1992; Grusius-Waetzel & 
Lesch 1998). Can magnetic reconnection be involved in 
the case of GRB 130925A? While this question remains 
open, certain aspects of it are already being addressed 
by ongoing studies (e.g., Sironi & Spitkovsky 2014) 

5.2. The origin of the radio emission and its 
connection to the X-rays 

The radio emission from most long GRBs is believed to 
originate from a forward shock in an external medium. 
If the radio emission of GRB 130925A is indeed origi¬ 
nating from a forward shock then the properties found 
in § 14.31 are the properties of the ISM (or GSM). We 
test whether the observed X-ray emission can arise from 
a simple forward shock afterglow model, using the best 
fit parameters we derived from the radio data. In this 
case, extrapolating the steep power-law (or mono ener¬ 
getic) synchrotron emission into the X-ray bandpass re¬ 
sults in a low flux, orders of magnitudes below the ob¬ 
served emission reported by Bellm et al. (2013), Piro et 
al. (2014), and Evans et al. (2014) (expected vF^, of 
^ 10-20 erg cm“^ s“^ compared to the observed one of 
^ 10“^^ ergcm“^s“^). 

Both Piro et al. and Evans et al. suggest that there is 
only a weak (or no) contribution from a forward shock to 
the observed X-ray emission. Instead they provide dif¬ 
ferent explanations for that emission, such as blackbody 
radiation or scattering of the prompt emission by dust. 
Both Piro et al. and Evans et al. conclude then that 
the GSM density must be very low, with n < 0.1 cm“^. 
In contrast, we find that, if the radio emission originates 
from a forward shock, then the lack of non-thermal X-ray 
emission is not due to low GSM density but rather due 
to the unusual steep energy distribution of the emitting 
electrons. In fact the electron density we find is higher 
by a factor of more than an order of magnitude than that 
of Piro et al. (2014) and Evans et al. (2014). 

Alternatively, it is possibile that the radio emission 
is the result of a reverse shock ploughing through the 
dense ejecta. In this scenario, the electron density we 
find is that of the expanding ejecta shell. This renders 
the comparison between the electron density we derive 
and the low density GSM environment found by Evans 
et al. and Piro et ah, irrelevant. In this scenario, the 
underlying weaker component that is observed in the ra¬ 
dio spectrum may originate from the forwards shock. In 
fact, a second component of F^, oc instead of a con¬ 
stant flat spectrum one, is also consistent with the data. 
Thus the existence of this second component makes the 
reverse shock scenario more plausible. However, if this 
second component is not from forward shock emission, it 
means that the afterglow from the forward shock is sup¬ 
pressed. As was already discussed by Evans et ah, one 
way to surpress forward shock emission is by invoking a 
very weak magnetic field in the forward shock, compared 
to the one in the shocked ejecta. Alternatively, this can 
also be achieved by reducing the kinetic energy to less 
than lO^^erg, an energy budget which is still in agree¬ 


ment with the minimum kinetic energy we find in § 14.31 
The reverse shock scenario, however, cannot explain the 
fact the we still see emission, above the secondary con¬ 
stant flux component, at 43 days. Any emission from a 
reverse shock is expected to decrease by at least 1 — 2 or¬ 
ders of magnitude compared to the emission observed at 
early times. Thus, our observations at late-times ham¬ 
pers the reverse shock scenario. 

The question, though, what is the origin of the X-ray 
emission, still remains. If we put the implications of the 
radio data aside, a forward shock as an explanation for 
the X-ray emission becomes unlikely also when compar¬ 
ing the emission at a wavelength of 2.2^m (Greiner et al. 
2014) to the X-ray flux, two days after discovery. The 
ratio between the two implies a relatively low spectral 
index between —0.3 to —0.6, while we expect a steeper 
spectrum due to the transition beyond the cooling fre¬ 
quency. In fact, the shockwave parameters derived by 
Piro et al. suggest a 2.2/im emission much higher than 
the observed one. 

It is possible, in some cases, to produce X-ray emission 
via the inverse-Gompton (IG) process. For example, a 
large enough reservoir of optical photons can be easily 
up-scattered to X-rays by electrons with a Lorentz factor 
of 7 e ~ 50. This possibility is especially intriguing since 
both the radio and the X-ray spectrum exhibit similar 
steep spectral slopes. Using the measurements of Greiner 
et al. as an estimate for the available optical and infrared 
photons, we find that the expected IC X-ray emission 
is vFi, ^ 5 X 10“^'^ erg cm“^ s“^. The observed X-ray 
emission, therefore, can not be accounted for by IG, since 
it is higher by almost two orders of magnitude than the 
expected IG emission. 

Another possibility is that the X-ray emission origi¬ 
nates from a different process and region in the progen¬ 
itor system, than the radio emission. As already men¬ 
tioned, Piro et al. suggest that most of the X-ray emis¬ 
sion is blackbody emission originating from a compact 
radius of 10^^ cm. This radius is much smaller than that 
of the radio emitting region. In contrast, in the Evans 
et al. model the X-ray echoing dust lies far away at pc 
scales. At the same time, one can argue that the fact that 
the emission at both radio and X-ray exhibit an unusual 
steep spectrum, may suggest a common origin. Addi¬ 
tional analysis of the observed X-ray emission is beyond 
the scope of this work. 

5.3. Implications for the nature of the progenitor 

As GRB 130925A belongs to the ultra-long GRB sub¬ 
class, we consider the two main progenitor scenarios sug¬ 
gested for this class, namely a TDE or an extreme col- 
lapsar. One of the best examples of a relativistic TDE 
is Swift! 1644-1-57 (Levan et al. 2011; Bloom et al. 2011; 
Burrows et al. 2011). Comparing the radio properties of 
Swift! 1644-1-57 (Zauderer et al. 2011; Berger et al. 2012) 
to GRB 130925A we find significant differences. First, 
Swift! 1644-1-57 was more radio luminous, by a factor of 
> 50 than GRB 130925A. Second, the radio spectrum 
peak of Swift! 1644-1-57 traversed below 10 GHz only at 
very late times, after 300 days. Furthermore, the exter¬ 
nal medium surrounding Swift! 1644-1-57, is much denser. 
Overall, the radio properties of Swift! 1644-1-57 do not 
resemble in any way those of GRB 130925A. This, of 
course, is not proof that GRB 130925A is not a TDE, 
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as Swift] 1644+57 is the first discovery of an onset of a 
TDE, but its properties may not be representative of the 
the whole population of relativistic TDEs. Beyond the 
comparison of the properties of GRB 130925A to the one 
of the Swift 16AA+57, the properties of the former do not 
meet the expectations of theoretical TDE models. Ac¬ 
cording to Giannios & Metzger (2011), the minimal rise 
time of the radio flux will be Gise ~ 30 — 200 days, de¬ 
pending on the black hole mass. In contrast, the peak 
radio emission of GRB 130925A is already declining be¬ 
tween 2 to 9 days, after discovery. 

The radio emission does not pose a clear challenge to 
the collapsar scenario. The radio can be explained as ei¬ 
ther a forward or a reverse shock in this scenario. How¬ 
ever, Piro et al. have suggested that the progenitor is a 
blue super giant (BSG) star with very little or no mass- 
loss. If the radio is originating from a forward shock, 
then a denser GSM than the one found by Piro et al. 
is required. Thus a progenitor with such low mass-loss 
as the one suggested by Piro et al. is not possible. On 
the other hand, the reverse shock case does not present 
a contradiction to the BSG progenitor scenario. 

In both of the above progenitor models, the spectral 
cut off in the radio, has to be explained, especially since 
it was not observed before in TDEs or normal long GRBs 
(collapsars). In fact, the unique radio spectrum that may 
require an alternative electron acceleration mechanism, 
may, in turn, point to a different progenitor system that 
has not been considered so far. For example, in § 15.11 
we speculated that magnetic reconnection may be the 
mechanism responsible for the acceleration of the radio 
emitting electrons. Singh et al. (2014) have studied the 
role of magnetic reconnection in accretion disk systems 
from micro-quasar scales to blazars. They found that fast 
reconnection originating from the central core, can play 
a part in the main energy output of these type of objects. 
At the same time, they found that this core mechanism 
cannot explain the energy output of GRBs and that the 
afterglow emission may arise from a more distant loca¬ 
tion, such as the jet. However, there is a suggestion by 
Giannios (2013) that magnetic reconnection may take 
place in jets as well. Gan it be that GRB 130925A is in 
fact not related to either a collapsar or a TDE event? 
Another piece of the puzzle is that the radio spectrum 
slowly decays to a flat spectrum in ^ 100 days. Can this 
be a hint that this event is actually related to an AGN 
activity? Assuming so, the minimum variability time 
of ~ Is in the prompt emission (Greiner et al. 2014), 
suggest a BH mass of < 10® Mq. This low mass disfa¬ 
vors the AGN scenario. Moreover, ~ 1.5 years after the 
GRB discovery, the seemingly constant flat radio emis¬ 
sion component has disappeared. 

6. SUMMARY 

We report here for the first time, the early ra¬ 
dio observation and detection of an ultra-long GRB, 
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